Gossamer metals 
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Laughlin's construction of exact gossamer ground states is applied to normal metals. We show 
that for each variational parameter < g < 1, the paramagnetic or ferromagnetic Gutzwiller wave 
function is the exact ground state of an extended Hubbard model with correlated hopping, with 
arbitrary particle density, non-interacting dispersion, and lattice dimensionality. The susceptibility 
and magnetization curves are obtained, showing that the Pauli susceptibility is enhanced by corre- 
lations. The elementary quasiparticle excitations are gapless, except for a half-filled band at g = 0, 
where a Mott transition from metal to insulator occurs. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.30.+h 



Progress in the understanding of many-body effects 
in strongly correlated electron systems, such as quan- 
tum magnets, narrow-band transition metal compounds, 
fractional quantum Hall systems, or high-temperature 
superconductors, has depended on a variety of theoret- 
ical tools. Important information about the electronic 
structure can often be obtained from ab initio calcula- 
tions, which are however less reliable if interactions be- 
tween electrons are dominant over their kinetic energy. 
On the other hand, the study of idealized model sys- 
tems, containing only the presumably relevant degrees 
of freedom, can provide insight into microscopic physical 
mechanisms. However, since such models are rarely ex- 
actly solvable, analytical and numerical calculations usu- 
ally involve approximations or extrapolations. In view of 
these limitations, support for proposed physical notions 
has occasionally come from an inverse strategy: start- 
ing from a correlated many-body wavefunction one con- 
structs a hopefully "reasonable" model Hamiltonian for 
which it is the exact ground state. Correlated quantum 
phases may then be classified according to their elemen- 
tary excitations or correlation functions. This approach 
has been useful in particular for the understanding of the 
fractional quantum Hall effect, spin-Peierls or Haldane- 
gap antiferromagnets, and quantum rotors 0. 

Recently, Laughlin Q developed a new approach to 
high-temperature superconductivity, viewing the insulat- 
ing state as a superconductor with very low superfiuid 
density. Pursuing the above strategy, he proposed that 
the ground-state wavefunction of such a "gossamer su- 
perconductor" is obtained from the BCS mean-field prod- 
uct state by applying the Gutzwiller correlation operator 
(0<5<1), 
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1 - (1 - g)Dt 
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where Di — n^n^ is the operator for double occupation 
at lattice site i, and constructed a corresponding model 
Hamiltonian. Elementary excitations [2L the transition 
from superconductor to Mott insulator [2j, magnetic in- 



stabilities 4| , and related mean- field Hamiltonians 
were also studied in this context. 

The purpose of this letter is the application of Laugh- 
lin's gossamer paradigm to normal metals, i.e., itinerant 
electrons on a lattice without broken (discrete) transla- 
tional symmetries. (In particular, antiferromagnetic or 
superconducting phases are excluded.) It is well known 
that a metallic system can be driven into an insulating 
state by strong electronic correlations. This type of tran- 
sition from metal to insulator, the Mott transition, oc- 
curs for example in transition metal oxides, and has been 
analyzed by a variety of theoretical methods @. These 
include the variational Gutzwiller wavefunction (GWF) 
obtained by acting with on an uncorrelated Fermi 
sea. In general the GWF describes a correlated metal, 
except for the insulating state with one immobile particle 
at each lattice site that results at g — for a half-filled 
band. When used as a variational wavefunction for the 
Hubbard model and evaluated within the Gutzwiller a; 
proximation Q, this Brinkman-Rice (BR) transition 
occurs at finite critical Hubbard interaction Uf R . While 
the Gutzwiller approximation becomes exact in the limit 
of infinite dimensions jSj, the BR transition is shifted to 
U® R = oo in finite dimensions [1(1 lllj . However, the 
reliability of these variational results is limited, as the 
true ground state of the Hubbard model in infinite di- 
mensions may behave rather differently ^3] i for example 
the number of doubly occupied sites in general does not 
vanish at the transition as in the BR scenario. Further- 
more, the analysis of elementary excitations is hampered 
by the fact that the true ground state is lower in energy, 
and on these grounds the GWF has been criticized as in- 
adequate for describing the Mott transition . Some of 
these difficulties are resolved for models with exact GWF 
ground states, which we now proceed to construct. 

Metallic gossamer ground state. In general a gossamer 
ground state is built as follows Q. Starting from an un- 
correlated product wave function \<j>) and operators b ka 
such that b 



= for all k and a, one applies an in- 
vertible many-body correlator K to obtain a correlated 
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wavefunction — K\<f>), and defines b ka = Kb ka K 1 . 
Then is an exact ground state of the hermitian Hamil- 

tonian H = J2kv E k<jb k(J b ka for arbitrary E ka > 0, since 
H > and H\ip) = 0. 

In the present context we use the Gutzwiller correlator 
(JTJ as in Refs. 0, 0, 0, IE which is invertible for g ^ 0, 
K(g)^ 1 — Klg^ 1 ), but start from a product state con- 
taining spin-up and spin-down fermions, characterized by 
the occupation numbers n kr7 (with n ka = or 1), 



n 



t,l°>- 



(2) 



ka « =1) 



This state is annihilated by the operators b ka = (1 — 
n ka)^ka n ka^ka- After some algebra, we can rewrite 
the Hamiltonian H as 



H = H t + Hh + H„ + H 



x 



Hy + const . 



(3) 



ffh = -/i^(n iT -n^), (4) 



(5) 
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#y = X! ^^snjaC^c^ , H v = U , (6) 

with the constant term depending only on g and the total 
particle density n = n-j- + hi, which is fixed; we mostly 
consider densities n < 1 since | -0 (0)} = otherwise. Here 
Tij a is the Fourier transform of E ka = (1 — 2n kcr )E ka , 
and the other parameters are given by 
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[(1 - g)T iia + (1 + ff)T <iCT 



(l + .g 2 )T yCT + (l- 5 2 )T, 



(7) 

(8) 

(9) 
(10) 



where Tij a is the Fourier transform of E ka , and L is the 
number of lattice sites. 

A model with arbitrary non-interacting dispersion e k 
can now be obtained as follows. For given band disper- 
sion E ka we construct the Fermi sea via n ka = 0(— E ka ) 
and let E ktT = \E ka \ > 0. Then we put E ka — e k - e Fa - 
and adjust the Fermi energies epa so that the starting 
wavefunction \<j>) is a Fermi sea with desired densities n a 
— j; ^2 k n ka . In the following we assume ^2 k e k = for 
convenience, hence = ^ fc e^rt^ < 0. 

For each < g < 1 the Gutzwiller wavefunction \i/j(g)) 
= K{g)\<f>) is the exact ground state of the extended Hub- 
bard Hamiltonian @. It contains the kinetic energy H t 



of a single band e k , which is independent of g, and an op- 
tional Zeeman term Hh, absent for h-\ — ny. For g < 1, 
H contains interactions that involve at most two sites: 
a repulsive on-site interaction Hjj and correlated hop- 
ping terms Hx and Hy, whose amplitudes are related 
by gY ija = (1 - g) 2 {Tija + X ija ). Note that X ija , Y ija , 
and U all diverge in the limit g — > 0. Similar interaction 
terms appear in models with superconducting gossamer 
ground states 0, 0, HJ , but here those states cannot 
be lower in energy than \ip(g)). Apart from g, the mag- 
netic field and the strength and range of the interactions 
depend on the chosen band dispersion e k and the densi- 
ties n„. To illustrate the behavior of the amplitude Ty CT 
appearing in (|§|)- (|TU)l we now discuss several examples. 

One- dimensional systems. The dispersion for a one- 
dimensional ring with nearest-neighbor hopping — t < 
is 6^ = — 2icos(fc). For the Fourier transform of E ka = 
kfe - crt I we find 



Tj±i,j a = t [2n a - 1 + i sin( 
At 



(11) 



Tj+rjcr = —r^ 7T W sin(7rn (T r) cos^n^) (12) 

n(r z — 1) 

+ cos(nn a r) sm(im a )] , |r| > 2, (13) 

which falls off algebraically at large distances. At half- 
filling (n a = 1/2) it is on the order of 1/r 2 and alternates 
in sign for even r, while vanishing for odd r. This long- 
range behavior of Tij a is rather generic. As another ex- 
ample we consider "1/r" hopping, Tj +r ,j a = it(—l) r /r 
with dispersion = tk, for which the corresponding 
Hubbard model was solved by Gebhard et al. (see Ref. 6 
for a review). We obtain 



T 



i-iyt 



1 — in(2n a — l)r — e 



-27rm CT rl 



(14) 



again with contributions proportional to 1 /r (absent for 
half- filling) and 1/r 2 . Similar power-law behavior is typ- 
ically found in dimensions D = 2, 3. 

Infinite- dimensional systems. Nearest- neighbor hop- 
ping t = l/\/2D on a hypercubic lattice with dispersion 
e k = —2tJ2 a=1 cosk a yields the density of states Phc(e) 
= exp(-e 2 /2)/v / 27r in the limit D -> 00 0,E3- In order 
to construct the corresponding amplitude Tij a further as- 
sumptions about its symmetry are necessary. Following 
Ref. [l8| we assume that it depends only on the "taxi-cab" 



D 



distance \ \R\ \ = J2 a 



ing Tija = T*J^2^) where r 



l tja -T,cr 

then obtain 



1 I R a I and use the appropriate seal- 

> 



1 / \ He r (e) 



T, 



2r+l,a - (1 - 2n ff )5 r , , 

He 2r (e Fo -) 



T 2r+2,a = 2phc(eRr) 



V(2r + 2)!' 



. We 

(15) 
(16) 
(17) 



3 




1 2 0.2 0.4 0.6 0.8 1 

h/t 1-g 



FIG. 1: Magnetization curves m(h) for a one-dimensional ring 
with nearest-neighbor hopping at half- filling. Inset: same for 
1/r hopping. 



FIG. 2: Expectation values of parts of H, for nearest- 

neighbor hopping in D = 1 (arrows) and D = oo at half- 
filling. Inset: double occupancy for 1/r hopping in D = 1, as 
compared to the pure Hubbard model. 



where He„(:r) are Hermite polynomials. At half-filling 
we find T 2 * ro . ~ r~ 5 / 4 , corresponding to an effective cor- 
related hopping range J2r T*a 2 °f order unity. For other 
densities of states p(e), in particular those with finite 
bandwidth, it is also possible construct a corresponding 
dispersion ek [l8j> an d then derive Ty- CT and Tij a in a 
similar fashion. 

Response to external magnetic field. Returning to the 
case of arbitrary dispersion and densities, we note that 
according to the equation of state J7J) the ground-state 
magnetization m — fit — fii is nonzero if an external 
magnetic field h is present. For the homogeneous suscep- 
tibility x we obtain 

xW' 1 = | L = iE 1 " (1 r"f K - ^ 

dm 4 ^ p(£Rr) 

In the limit of zero field this reduces to x(0) = Xo/[l — 
(1— g 2 )n/2]. As expected the system behaves like a corre- 
lated paramagnet, i.e., the interactions enhance the Pauli 
susceptibility xo = ^p( £ f) of the uncorrelated system. 
However, it should be kept in mind that the interac- 
tion parameters (|% ]) - (|10[1 do not remain constant when 
the parameters h or m are varied. Fig. ^ shows the 
magnetization as a function of magnetic field for a one- 
dimensional ring at half-filling. Interestingly, for nearest- 
neighbor hopping the upward curvature of these magne- 
tization curves is very similar to Bethe-ansatz results for 
the pure Hubbard model 0,H3> where a metal-insulator 
transition occurs at U c = + . By contrast, for 1/r hop- 
ping the magnetization curves are strictly linear, m = 
x(0)/i, due to the constant density of states. 

Metal-insulator transition. For an unpolarized half- 
filled band (n = 1, e~p a = 0), the ground-state wave func- 
tion \ij}(g)} describes a metal for g > and an insulator 
for g = 0. In the insulating state there are no doubly 



occupied sites, the discontinuity of at the Fermi sur- 
face vanishes, and the kinetic energy (H t ) is zero. This 
Mott metal-insulator transition in the ground state of H 
occurs at infinite interactions (|5|)- (|10fl . in contrast to the 
variational BR transition, or numerical results for the 
Hubbard model in infinite dimensions [12(. 

Nevertheless we may, somewhat artificially, shift the 
transition to finite interactions as follows. Clearly \ip(g)) 
remains the ground state when we multiply if by a pos- 
itive ^-dependent factor, although qualitatively different 
Hamiltonians may then result in the limit g — ► 0. For ex- 
ample, for the Hamiltonian = gH the X term has a 
finite limit, while = g 2 H yields a vanishing X term 
and finite Y and U terms; in both cases the quadratic 
kinetic energy vanishes at g = 0. In particular we may 
conclude that for any dispersion the Hamiltonian 

H' = Y l] il isn ]S c^ a c- a + U'Y^ n»t«il > (19) 

where is the Fourier transform of efc(l — n^ a ), has 
the exact ground state \ip(g = 0)) at half-filling if U' > 
U' c , with critical interaction U' c — — J2 a e 0o- = l £ o|- In- 
terestingly, the uncorrelated kinetic energy also sets the 
energy scale of the BR transition in the Gutzwiller ap- 
proximation H, where Uf K — 8 jeo j • Although H' is not 
a standard Hubbard Hamiltonian, it nonetheless appears 
to be the simplest model with a BR-type transition to an 
exact insulating ground state at finite Hubbard interac- 
tion. 

Ground-state expectation values. The separate expec- 
tation values of the kinetic energy H t and the interac- 
tion terms Hx, Hy, and Hu can be calculated from the 
quantities n ka = (c^ CT c fc(T ), d = ^.(".t^ll- x ij<r = 
("iaC^c^J, and y i]a = {h^n^c^c^) . Using the meth- 
ods of |3 0] it can be shown that for the GWF the 
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Fourier transforms of the latter are given by 



Vka 



1 - nka 



1 



1 



~9 

nka 



(i ~ n °ka)g n k* 

(i - n L).9 2 "^ 

fl - a) 2 



d, (20) 



-d,(21) 



i.e., only the momentum-space distribution rik a and the 
double occupancy d are needed. They may be evaluated 
numerically by Monte-Carlo methods, but are also avail- 
able in closed form under certain circumstances. For one- 
dimensional systems with symmetric Fermi sea (n kcr = 
n° , ) both quantities have been calculated analytically 
|l4l Il5| . In dimension D = 2, 3 high-order perturbative 
methods can be used [llj . The Gutzwiller approximation 
, with piecewise constant momentum distribution n\~ a , 
is recovered in infinite dimensions 0- 

The expectation values of the various parts of H are 
shown in Fig. [21 for nearest-neighbor hopping in D = 1 
and D = oo at half- filling. We note that (Hx) approaches 
a constant for g — > 0, while (Hy) and (Hu) diverge. This 
behavior occurs for all D, since d ~ g 2 ln(l/g) in one di- 
mension d = o(g) in all finite dimensions [ToL Hl|. 
and d ~ <? in infinite dimensions. We may thus conclude 
that the penalty that Hjj imposes on double occupancies 
is compensated by assisted hopping due to the nonstan- 
dard three-body interaction Hy. 

The effect of correlated hopping is also apparent when 
comparing to the pure Hubbard ring with 1/r hopping, 
which features a metal-insulator transition at IL- = 2nt 
with continuous nonzero double occupancy d 6]. For 
comparison with previous studies of variational wave- 
functions in the vicinity of this transition [Hi. Il7j . d vs. 
U is shown in the inset of Fig. [21 The results for both 
models with 1/r hopping agree for weak interactions, but 
the energy gain from correlated hopping leads to a larger 
number of doubly occupied sites for strong coupling in 
the model ©, as expected. 

Quasiparticle excitations. The known ground state of 
H suggests that it might also be possible to calculate 
dynamical properties of the model, such as the spectral 
function. Unfortunately the construction of exact excited 
states is not straightforward, be it with one added or 
removed particle, or with charge or spin excitations. We 
therefore proceed by considering the variational states 

"El 



Kc + ka \ 



Kc 



ka\ 



whose mean energy is 

(ka\H\ka) 



El 



_ (V#) 

(ka\ka) (ka\ka) 



■E k0 



(22) 



(23) 



where the commutator relations [b k , bt, ,] — Skk'daa' 
and [b kcr , b k , a ,\ = were used. The states \ka) are mutu- 
ally orthogonal and their energy is thus an upper bound 




FIG. 3: Quasiparticle excitations in a one-dimensional ring 
with nearest-neighbor hopping t > 0. 



to the quasiparticle energy for momentum k and spin a. 



The variational energy to add a particle (i.e. 



E + = l £ fc ~ £ Rt| 

ha 1 - (1 + g)[(l - g)n s + (1 + g)n 



while for the removal of a particle (with n£ CT = 1) 



E. 



9 r\e h - e F(T | 



kcr 



(1 + - 9)n a + (1 + g)n ka ) -l-2g 



0) is 
(24) 



(25) 



Clearly the quasiparticle excitations are gapless, since 
E k(T — > close to the Fermi surface. Fig. shows these 
energies for one-dimensional nearest-neighbor hopping at 
half-filling. 

Conclusion. We have constructed and characterized 
a new class of itinerant electron models for which the 
metallic Gutzwiller wavefunction is an exact ground 
state, due to the interplay of Hubbard interaction and 
correlated hopping. For a half-filled band a Mott metal- 
insulator transition similar to the Brinkman-Rice sce- 
nario occurs, illustrating Mott's original idea of a quan- 
tum phase transition entirely due to charge correlations 
without magnetic ordering. Further study of the elemen- 
tary excitations in these models should be fruitful. 

This work was supported in part by the DFG via 
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